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Abstract 

We present a study of the Drude weight D(T) of the spin-1/2 XX Z chain in the gapless 
regime. The thermodynamic Bethe ansatz (TBA) is applied in two different ways. In the 
first application we employ the particle basis of magnons and their bound states. In this 
case we rederive and considerably extend earlier work in the literature. However, in the 
course of our investigation we find arguments that cast doubt on the applicability of the 
TBA in this case. In a second application by use of the spinon and anti-spinon particle 
basis we obtain completely different results. Only for anisotropy parameter A close to 
we find that D(T) is a monotonously decaying function of temperature. For A close to 1 
the behaviour is entirely different showing a finite temperature maximum. Also for the 
isotropic antiferromagnetic chain (A = 1) the results for D(T) are finite for T — as well 
as for T > with an infinite positive slope at T = 0. 

1 Introduction 

In this work we are studying transport properties of the spin-1/2 Heisenberg chain with 
longitudinal anisotropy (XXZ chain). Depending on the representation of the system the 
quantity of interest is the spin conductivity of a quantum spin system with anisotropic spin 
exchange or the electrical conductivity of a system of spinless fermions with density-density 
interaction. Both representations are related by a Jordan- Wigner transformation. The re- 
cent interest in these quantities has several reasons. We want to mention only two of these. 
First, in low-dimensional systems the question has been raised whether spin diffusion exists 
or not and the role of integrability for anomalous transport properties was discussed. (For a 
review see pfl.) Second, the Drude peak at zero frequency in the dynamical conductivity is 
in principle accessible to analytical studies. 

Let us consider a Hamiltonian of the form 

H = H + V = -tJ2 (e ieA ^ct +lCj + e- ieA ^c]c j+ i) + H int (1) 

3 

representing a one-dimensional system of length L with periodic boundary conditions subject 
to a vector potential A. In our case Hint describes density-density interactions and hence 
does not depend on the vector potential. 

Within Kubo theory [2], i.e. linear response for (QJ) in A, the dynamical conductivity is 
obtained in terms of current-current correlation functions. The Drude weight D is the zero 
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frequency distribution of the dynamical conductivity a = D5(lu) + <r rcg and has a spectral 
representation in terms of eigenstates and eigenvalues of the system with vector potential 
A = 



2L 



n\j\m 



(2) 



where j is the current operator, (Tj is the thermal expectation value of the kinetic energy, 

and p n = e~ tn / T /Z is the Boltzmann weight for the eigenstate \n) of the Hamiltonian with 
energy e n . 

On the other hand, a static magnetic flux leads to a site and time independent vector po- 
tential A x (j,t) = <fi/e with characteristic dependence of the eigenvalues on <j). Denoting the 
Hamiltonian Q by H(4>) we find in second order pertubation theory in <j) 



e n (<j>) = (n\H(0)\n) -</)(n\j\n 
Comparing second order terms we see that 



n\j\m 



n\T\n 



(3) 



D 



1 ^ 1 



1 d 2 e r , 



(4) 



This is the generalization of Kohn's result see also [UEj to finite temperature In this 
way the Drude weight D is connected to the twist <f> in the boundary conditions caused by 
the applied external field. 

The expression (J3J is very interesting as it allows for the calculation of the Drude weight just 
from the eigenvalues of the Hamiltonian without the knowledge of matrix elements. This is 
essential for analytic calculations of integrable systems rendering the task feasable although 
the remaining work is still formidable. In [Jj an extension of the traditional Thermodynamic 
Bethe Ansatz (TBA) to cover the mean curvature of energy levels was presented. This method 
was applied to the spin-1/2 XX Z chain in |5J. The numerical evaluation of the Drude weight 
showed a curious behaviour especially for the isotropic chain suggesting that D(T > 0) = 0. 

Our own interest in the finite temperature Drude weight was triggered by the peculiar findings 
of |Hj. The treatment of the (strictly) isotropic Heisenberg chain within the TBA approach is 
rather challenging as it involves infinitely many functions to be solved from infinitely many 
non-linear integral equations. Hence, we suspected that the findings in |Sj might be based 
on inappropriate numerical treatments of these integral equations. This however, is not the 
case as will become clear in section 3. 



In Sec. 2 we present some necessary elements of TBA and the generalization for calculating 
the Drude weight along reference 0. In Sec. 3 we describe the application of the TBA 
approach to the Heisenberg chain on the particle basis of magnons and their bound states 
(strings) somewhat following the treatment of 8 . Our main technical achievement in this 
section is that we manage to reduce the many integral equations to just two. Among other 
things, this allows for much simplified numerical calculations in comparison to jHj, however 
confirming the obtained data and the steep drop of the Drude weight as discussed above. 
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In section 4, we apply the extended TBA method on the basis of spinons and anti-spinons. 
The resulting equations are studied analytically and numerically yielding results that are 
completely different from those of section 3 and ref . jS] . Most of the discussion of our results 
is presented in section 5. Some more technical material related to the derivation presented 
in section 4 can be found in the appendix. 



2 TBA formalism for free energy and Drude weight 

We consider a system of particles with bare energy 6 a (x) and momentum Pa{x) parametrized 
by the spectral parameter x. The index a labels the species of particles. As we do not use 
in this work the notion of the dressed energy function we take the liberty of dropping the 
commonly used upper index (e a (x) instead of e« (x)). 

Furthermore we have diagonal scattering of two particles of species a and (3 with scattering 
phase Q a p{x a — xp). The quantization condition for an eigenstate characterized by a set of 
particles with rapidities x a k for twisted boundary condition with angle (/> reads 

Lp a (x ak ) + ^ ®a(3{x a k ~ %/3l) = 2lll ak + (j) a , (5) 

pi 

where <p a is a multiple of the applied twist angle <p a = 4> ■ n a with some (integer) number 
n a . (Note that Bethe ansatz equations usually take such a form.) We introduce the counting 
function Z a 

The set of solutions to ® for x ak in an interval of width Ax a comprises actual rapidities as 
well as holes with density functions p a and p^ the sum of which is related to the counting 
function Z a 

^particles + Wholes = (p a + p^)Ax = Z a (x + Ax) — Z a (x), (7) 

or explicitly 

Pa (x) + P h a { X ) = ^Pa(x) + ^Y,Y1 @ '^ X ~ X Pti- ( 8 ) 

B I 

The summation over rapidities on the r.h.s. can be written in the thermodynamic limit in 
terms of integrals over pp 

1 1 , 

(l + Va)pa(x) = — p' a {x) + — J2x a p*pf3, >i a B := &' a(3 , T} a \= y . (9) 

This set of integral equations is equivalent to the Bethe ansatz equations (JSJ). The state 
representing the macrostate for finite temperature T is obtained from the minimization of 
the free energy functional 

/ = e-Ts 

= y~*A e a p a dx-T ( [p a + p^] \n(p a + p h a ) - p a lnp a - p^lnp^j dx). 

(10) 
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This results into the non- linear integral equations (thermodynamic Bethe ansatz equations) 

lwq a (x) = f3e a {x) - — ^ x a p * In (l + n^ 1 ) , (11) 
77 P 

where we have used the symmetry property K a p(x) := xp a (—x). By use of this set of 
equations we can simplify the expression for the free energy ()1U|) yielding 

1 f°° 

The equations (|10|) - (jl2|) are valid for finite magnetic field h if the term e a is replaced by 
e a — hn a where n a is the same number that occurred below We do not give the details 
of these calculations. The interested reader is referred to the book 

Here we are more concerned with the finite size analysis of the rapidities for which we closely 
follow the treatment of [7j. For x a j we make the ansatz 

JX) „( 2 ) JX)( oo\ „( 2 )/ oo\ 

(1 2*1 (1 2) 

with finite size coefficients g a - taking the form of smooth functions g a ' (x) in the thermo- 
dynamic limit. Inserting this into the counting function leads to an expansion 



P I 



. g pi 9c 



P l 



+ ^£E°V*-*5>(i) 
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" iEE^(x-x^)|. (14) 
We identify the O(l), 0(1/L), 0(1/L 2 ) contributions to the counting function 

7 (2) 

Z B (x) = 2S°(x) + -^(x) + -g-(x) (15) 



4 1} (x) = -^E '^*^)-^ (17) 

p 

= ^E(l @ U(af P p)-e' aP *(9f ) Pp)), (is) 
p v y 

where we have replaced the summations over x^ by integrals involving the density functions. 
For studying the quantization condition Q we have to expand terms of Z like 

/ 9 {1) 9 (2) \ 

Z a (x aj ) = zJx^ + y -f + y -^\. (19) 
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leading to 



Z a (x a j) 



+ 



L 2 



u aj >i>aj 
1 

I 



+ 



L 2 



(i) 2 

aj n>aj 



+ Z™ { x a°j)9aj 



2irl n 



This equation imposes for the corrections 

Zj»{x%) + ZZ'(x%)gW=0 
Z< 2 >(*~.) + zW(x%)g<» + \zf(x^)£f + Z?(x%)g$j = 0, 



reading in the thermodynamic limit 

(1 + Va) ■ {9$ Pa) = 27T + ^ * ^fpp) 

(1 + Va) ■ (g ( a ] Pa) = l' a + ^ ^ X «£ * (ffl Pfi) 



la = \(Pa + Pa)g ( a )2 ~ ^ Yl ^ * ^ Pp) 



Focusing again on the Drude weight we are rather interested in the derivatives with : 
to 4> which we denote by dots 

(1 +Va)- (9^ Pa) = ■^ + ^J2 >< aP*(.9 i f } ) pp), 

(1 + Va) ■ (9a 2) Pa) = l'a + ^zY,^aP*{gfpp), 

Z7T P 

la = {Pa + Pa)g ( a )2 ~ 7TL Yl K ^ * ^ Pp)- 

Z7T P 

Next we note the second order term in <fi of the energy function 

E2 = Y,j (l e *"9a 1)2 Pa + ea'g a 2) Pa^J dx, 

the second derivative of this is 

/oo 
Ua"g a 1)2 pa + e a 'g^p a )dx. 
-00 ^ ' 



At this point we want to comment on some potentially fatal problem residing in the last 
expressions. By inspection of explicit examples we see that separate integrals of the individual 
summands occurring in the integrands of l|24jl and (|25|) diverge! 
Ignoring this problem we may reformulate the last expression 

1 / /'OO /*oo \ 

D = 2^(j ^"(9 { a )2 pa)dx + j e a 'gW Pa dx) 

= w^L dx > (26) 

where in the last line we have used among other things an identity following from the "dressed 
function" formalism applied to (after taking the derivative with respect to the the spec- 
tral parameter) and (|24j) . 

Expression (|2fij) was derived in [Jj for the study of the Hubbard model. Here we like to 
simplify this expression by showing that all functions appearing in the integrand of ()26|) are 
simply related to rj a . For this we note the relation 

ea = J^ P ' a , (27) 

7 

and the validity of (fTT|) for finite field h provided e a is replaced by e a — hn a . This yields 

9^ Pa = i^^l^l + r?- 1 ), (28) 

just because both sides of each equation satisfy the same integral equation. Putting together 
422) and (JUD we find 

D = Jsm7 v r y^^jj^fag dx 

4tt/3 7 ^J-oo (1 + %)(! + r/a 1 )^ hi % 



This expression is amazingly symmetric as the identity 

(^Inr^flnr?) 2 In In ± 



(30) 

(1 + t/)(1 + r/- 1 )^ lnr/ (1 + ^(l + r,)%, In ± 

may be interpreted (up to the sign change) as invariance under a "particle hole transforma- 
tion" (i] = p h /p <-> 1/rj = p/p h )- Formulas like (j29j) with the additional minus sign also 
appeared in [S]. 



3 TBA and Fusion Hierarchy for the Heisenberg chain 

In general the number of particles entering the TBA formalism as sketched in the foregoing 
section is infinite. In the case of the spin-1/2 Heisenberg chain we have to deal with the single 
magnon and its many bound states (strings). For the special case of 

H = jJ2[SfSf +1 + S*Sy +1 + co8( 7 )5? S? +1 ] -hJ^Sf (31) 

i i 
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with anisotropy parameter 7 = tt/u with integer v the number of bound states is finite and 
the TBA equations close for a finite set of 77-functions. This is the case studied in |Hj with 
the main equation being (|26j) for which the functions rj etc. are calculated from (jllj) etc. 

We want to calculate the functions appearing in (|29[) in a different manner by employing 
an alternative approach to the thermodynamics of quantum systems making use of a lattice 
path integral formulation |1U1 lllj . The quantum system at finite temperature is mapped to 
a classical model. In the case of the Heisenberg chain we are led to the study of a staggered 
six-vertex model on the square lattice. The size of this lattice is L x N where L is the 
length of the quantum chain and N is the Trotter number. The partition function of this 
model is calculated within a transfer matrix approach. For these calculations it turns out 
that the quantum transfer matrix (QTM), i.e. the column-to-column transfer matrix, is most 
appropriate as it shows a spectral gap between the largest and next-largest eigenvalues even 
in the limit L, N — > 00. 

There are mainly two different ways of analysing the spectrum of the QTM: the Bethe ansatz 
analysis and the fusion hierarchy. Here we make use of both approaches. The fusion hierarchy 
was studied carefully in which we are going to utilize extensively in the first part of this 
section. In a hierarchy of functions Y a (x), a = 1,2,..., is derived from the basic object 
of the QTM. These functions satisfy the functional equations 

Y a (x + i)Y a (x - 1) = (1 + r a -i(a:))(l + Y a+1 {x)), (32) 

relating each function Y a to two other functions, one with lower index, the other with higher 
index. For anisotropy 7 = ir/v with integer v the functional equations close at a finite level. 
(In ^2j v is denoted by po; otherwise we follow their notation closely.) The closure happens 
as 1 + Y u -i factorizes 

1 + (x) = (l + ^ hu/2 K{x)} (l + e"^/ 2 K(x)) (33) 
with a function K (x) satisfying 

K{x + x)K{x-i) = l + Y v - 2 {x). (34) 
These equations can be put into a more canonical form by the definition 

Va(x) '■= Y a (x), for 1 < a < v — 2 
^_i(x) := e + ^' 2 K{x), 

rj v {x) := e +phlj/2 /K(x). (35) 

These i/-many functions satisfy the functional equations 

rj a (x + i)rf a (x - 1) = (1 + 77a-i(x))(l + r) a+1 (x)), for 1 < a < v - 3 
l)v-l{x + i)Vv-2 (x ~ 1) = (1 + r/„_ 3 (x))(l +rj u -i(x))(l +rj~ 1 (x)) 
e-^V-i(z + i)^-i(x-i) = e + ^ hu rj- 1 (x + i)7 1 ~ 1 (x-i) = l + 7 lu . 2 (x). (36) 

These equations can be brought into the form of z^-many non-linear integral equations iden- 
tical to the TBA equations of 9 that are the starting point of jSj- 

We want to present a method of computation of the relevant functions avoiding an explicit 
solution of the TBA equations. To this end we observe two important properties of the func- 
tions rj a : 
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(i) All objects Y a and K are even functions of the magnetic field h. Hence the first order 
derivatives of rj a with respect to h evaluated at h = yield zero for all a except a = v — 1, 
v where we have 

d 0h log 77^-1 = dph log i] u = v/2. (37) 

(ii) For vanishing magnetic field we have a remarkably simple identity 



1 + T}„-i(x) = 1 + r] u 1 (x) 

Q{x + i{v-l)Q{x-i(v -I) 

[sinh (^(x + i(u + v — 1))) sinh (^(x — \{u + v — 1)))] 



W/2 



(38) 



involving the function Q that can be computed without solving the z^-many TBA equations. 
The alternative computation is done on the basis of a set of non-linear integral equations for 
two functions a and a 



logo(x) = -/3e + (x) + I k(x - y) log(l + a(y))dy 

k(x -y + 2i) log(l + a(y))dy 
loga(x) = -/3e-(x) + I k(x - y) log(l + a(y))dy 



-oc 
oo 



-oo 

^DO 



J — OO 

/oo 
/ t (x-y-2i)log(l + a(y))dy, (39) 
-oo 



where the driving terms e± and integration kernels k are explicitly given by 

7TSin7 7T 1 

e± x = J— e (x ±— -/i, e (x = — — , 40 

2 7 2(7r — 7j cosh|x 

1 I'" sinh(^-2)fc 
= — / 2— ? ^—e tkx dk . (41) 

27T ' 



oo 2 cosh /c sinh ( - — 1 ) fc 



In the course of the derivation of these two non-linear integral equations, explicit expressions 
for the function logQ(x) occur. By use of these results we find for (|3*5|l straightforwardly 

1 + 7]„-i{x) = 1 + r/~ x (x) 

= zv-exp^y [w(x-y-i)log(l + a(y))-o;(x-y + i)log(l + a(y})]d^ (42) 
where the function w(...) is given by 



u( x ) = ^ / y~ ^e ikx dk . (43) 

-oo 2sinh [ z - 1 ] k 



2vr 
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The /c-integral may be regularized by either choosing a path along the real axis avoiding k = 
in the upper or in the lower half plane: the term in brackets on the r.h.s. of (|42|) containing 
convolution terms of uj with a and a does not depend on the choice of regularization as the 
asymptotics of a and a are identical. 

The last equations are the main result of this section. The r.h.s. of ()29|) can be calculated as 
the terms with 1 < a < v — 2 are zero and those with a = v — 1, v are obtained from (|37|) 
and (02). Most importantly, in this formulation the number u, that used to be the level of 
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closure of the TBA equations, only enters via the anisotropy parameter 7. Hence, we may 
regard the above formulation as an analytic continuation to all anisotropies < 7 < ir. 
Finally, we like to mention that the largest eigenvalue may be calculated from 



where an irrelevant constant, i.e. independent of T and h has been ignored. 
3.1 Numerical results 

Here we present the results of numerical treatments of the equations (|39|) - (|42j) derived above. 
For any anisotropy < A < 1 we obtain Drude weights D(T) that are monotonously decaying 
for increasing temperature T, see Fig^ In the zero temperature limit the analytically known 
values -Do are reproduced (the formula is given in (|H4[) below). At high temperatures the 
data D(T) follow the behaviour derived analytically in (|48j) . 

Most striking is the behaviour of D(T) for parameters A close to 1. At T = the slope of 
D(T) gets steeper the closer A approaches 1. For A = 1 the Drude weight shows singular 
behaviour with a finite value of D(T = 0) and a drop to D(T) = for all T > 0. This is in 
perfect agreement with the results obtained in [Sj by a completely different and much more 
elaborate treatment of the TBA equations comprising v many non-linear integral equations 
(with v — > 00 for A — > 1). 

On the technical side we like to note that our treatment is not restricted to anisotropy param- 
eters A = cos ^ with integer v. Also, we conclude from the agreement of the two different 
numerical treatments of the TBA equations (here and in |Sj) that both results are either 
correct, or both are wrong in which case the failure must arise from assumptions underlying 
the analytical derivation of (|26|). We will return to this discussion at the end of the next 
subsection. 



3.2 High temperature asymptotics 

The analytical solution to the non-linear integral equations (NLIE) (jSHJ in first order in 
(3 = 1/T is obtained by writing the auxiliary functions 



These expressions are inserted into the NLIE and log(l + a),log(l + a) are expanded up to 
first order in (3 leading to the linear integral equations 




(44) 



log a(x) 
log a(x) 



0e{x) + 0(p 2 ) 
pe{x) + 0{(3 2 ). 



(45) 



/3e(x) 



00 



00 




k (x - y + 2i) e(y) dy 



—00 



—00 



-pe(x) 



/9e_(x) 



CO 



00 




—00 



—00 



(46) 
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T/J 

Figure 1: The Drude weight in the temperature range T/J = 0, 1 for different anisotropy 
parameters A = 0, 0.1, ...,0.9 as obtained in the TBA approach. For A = 1 the corresponding 
data consist of a finite value of the Drude weight for T = and zero for all T > 0. 



This system of linear integral equations is solved in Fourier space yielding explicit expressions 
for the Fourier transforms of e and e. Transforming back yields the result 

loga(x) = — sin 7 Bi ( - ; 5— )— Bh. 

y ' ' " \ \ _ e — 7x— 71 y — g— 7^+71 1 — g— i x ~ 371 / 

loga(x) = — shx^Bii- h ; — ) + Bh 

v ' ' r V \ _ g— 7X— 71 Y — g— 7^+071 I _ g— 7X+71 / ' 

(47) 

This has to be inserted into the expressions (|42|) and (|29j) for the Drude weight D. The 
integral was solved by means of Mathematica yielding the analytical high temperature result 

D ^ c(A) = J *2dL^l, m 

for 7 = jj, v = 2, 3, .... This was compared to the numerical data for D(T)T obtained by 
numerical iteration of the NLIEs for several anisotropy parameters in the high temperature 
limit. We observed very good agreement in this regime with high accuracy. (Prior to publi- 
cation of our work we communicated the analytical results with the authors of who also 
found very good agreement with the numerical data of jSj.) 

A fundamental problem appears with (|48|) if we recall the unitary transformation U acting 
on every second lattice site j by a rotation in spin space by an angle ir around the z-axis and 
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thus transforming SJ, Sj, 5| into —Sj,—Sj,Sj. Actually, in fermionic operator language this 
is equivalent to changing c^,Cj to —c^—Cj on every second site without affecting the anti- 
commutation rules. Hence we have UH( J, A)J7 _1 = H(—J, —A) and from (JIJ we conclude 

D(T, J, A) = D(T, -J, -A) = -D(-T, J, -A), (49) 

where we also used the fact that U does not change the dependence of the energy levels on 
the twist angle <p. 

Equation (|49|) relates the Drude weight at positive temperature T and interaction parameters 
J, A to the Drude weight at negative temperature — T and interaction parameters J, —A. In 
general we do not learn much from this with respect to the properties of a system at "phys- 
ical" temperatures, e.g. for T — > 0+ equation (|i9"|) relates the Drude weight of the ground 
state of the system with J, A to the Drude weight of the highest energy state of the system 
with J, —A. Quite differently, for high temperatures the asymptotical behaviour (|48|) is 
valid for large positive and negative temperatures T! In fact, a high temperature series in 
(3 (= 1/T) is meaningful and has to have some finite convergence radius such that sufficiently 
small, but otherwise arbitrary (complex) arguments (5 are allowed. Hence we may apply (|48|) 
in the form D ~ C(J, A)/T for positive and negative temperatures and from (|49|) we obtain 
C(J, A) = C(J, — A) or C(7) = C(ir — 7) a relation certainly not satisfied by (|48|) . There 
are several possible reasons for this failure. It may simply be that (|48j) is valid only for the 
discrete values of 7 = ^ with v = 2,3,.... In other words, the analytical continuation is 
simply not allowed (at least not in the "most natural" way). We do not want to discuss this 
point any further. 

Another reason may be due to assumptions usually taken for granted in TBA calculations, 
but failing in the derivation of (|26j) . We already pointed out that there are convergence 
problems in the transformation of (|25j) to (|2fi[) . Yet another problem is due to the fact that 
(|26j) has been applied to the Heisenberg chain by dealing with magnons and their bound states 
("string") as elementary, i.e. stable particles. This is -as we nowadays know- not strictly the 
case. The perfect string picture may be strongly violated, especially if the string is located 
at large spectral parameters. This explanation is actually corroborated by numerical studies 

4 TBA based on spinons and antispinons 

The above set of NLIEs (|39[) has been derived in the QTM setting which is mathe- 
matically quite different from the TBA. In this sense the QTM approach is an independent 
method based on algebraic and analytical reasoning rather than combinatorial arguments in 
the case of the TBA approach. However, the structure of the NLIEs 1)39(1 is similar to that 
of TBA. In fact, it has been argued Jl] that (|39|) can be viewed as the TBA equations of 
spinons and anti-spinons with bare energies and scattering data given in (J41[) (with 771 := 1/a 
and rj2 '■= l/o). Such a description of the thermodynamics of a system based on its ex- 
act low energy excitations would usually be expected to be possible, but restricted to the 
low temperature regime. Surprisingly, for the thermodynamical potential of the Heisenberg 
chain it is quantitatively correct for all temperatures and fields! Note also that ()39|) is valid 
for all 7, not just for the discrete set for which the number of bound states of magnons is finite. 

The idea of our work in this section is to utilize the spinon and anti-spinon particle basis for 
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describing the spectrum of the model. The advantage in comparison to the basis of magnons 
and their bound states is obvious: from the beginning there will be just two coupled NLIEs 
to be solved and only two terms in the sum of 1)26(1 in contrast to in general infinitly many 
for generic values of A. Also we no longer deal with disintegrating or deformating strings. 
Of course, this approach is phenomenological and for the computation of the Drude weight 
it can not be expected to be as successful as for the free energy. With respect to the justifi- 
cation of our approach and its limitation to low temperatures some arguments are given in 
the appendix. 

The reader short of time may directly go to equation (|56|) obtained from 1)29)) for just two 
particles: spinon and anti-spinon with r}\ := 1/a and r]2 := 1/a. In the following treatment 
we want to address certain subtle, but important issues like the response of spinons and anti- 
spinons to twist angles <j> and magnetic fields h. In section 2, the two types of response were 
assumed for any particle species a to be governed by just one number n a which is usually 
an integer. From (|41|) we see that the Zeeman energy for spinons involves a number ± 2(71-^-7) 
that in general is non-integer. From JH] we may infer that also the response to the twist 
angle (ft is governed by the same number. Here however, we want to treat the problem in 
a different way by a mapping to a system where the numbers n a are just 1 as usual. This 
system will involve just one closed integration contour thereby also avoiding the divergence 
problem mentioned in the final paragraph of the preceding section. 

To achieve the outlined goals we reformulate the equations 1)39)1 from which we see that the 
functions a(x) and a(x) can be analytically continued into the complex plane where they turn 
into each other due to the identity a(x) = l/a(x — 2i). 

Substituting log(l + a) = log(l + 1/a) + log a in (|3"5)l and resolving for log a yields 
log 77+ (x) = 0e(x) - I x{x-y)\og(l + r]~ 1 {y))dy 



j —00 

/oo 
■ 
-00 



x{x-y + 2i)log(l+r l I 1 {y))dy, (50) 

where 77+ (x) := a{x) and r]-(x) := l/a(x). The function e(x) happens to be the energy of a 
single magnon parametrized by the spectral parameter x 

e(x) = —J -p'(x) — h, 

7 

V ( X ) - n sinh^QE-i)) 

p{x) - ^og sinh( 2 (x + l)) ' i 51 ) 
and p(x) is the momentum. Furthermore, k{x) is given by 

«(,) = £*(„, e W; =-n„ g ^||^, (52) 

where 0(x — y) is the scattering phase of two magnons with spectral parameters x and y. 
Due to the relation r/_(x + 2i) = r] + {x) we can immediately write down a second equation 
from 1)50)1 yielding a complete set of equations for r/±. 

These equations are quite different from the standard TBA equations for the spin- 1/2 Heisen- 
berg model since those are based on up to infinitely many bound states (strings). However, 
the above equations may be obtained within the TBA method if the following two modifica- 
tions with respect to the standard approach are applied 
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• ignore bound states, 



• consider the single magnon with energy e(x) and momentum p(x) parametrized by the 
spectral parameter x on the real axis Im x = ("C+") as well as on the axis Im x = —2 
("C_"). 

A posteriori we find that this prescription renders the single magnon a complete particle 
basis for the thermodynamics of the system. In other words, if we take the magnon with bare 
momentum p and energy e on the axes C± the TBA method of the previous section yields 
exactly 1)50(1 , The free energy then reads 

1 r°° 

= Y £ / p '^ x) ln(1 + dx - (53) 

a=± J ~°° 

where p+(x) = p(x), P—(x) = —p(x — 2i) (both functions are monotonously increasing with 
increasing x) and r] a are determined from (|50|) or 

lwn a (x) = 0e a {x) - — ^ x a/3 * In (l + rip' 1 ) , (54) 

with e + (x) = e(x), e_(x) = e(x — 2i) and x ++ (z) = x (z) = x(z), X-\ (z) = —x(z + 2i), 

x |_(z) = — x{z — 2i). 

Finally, we note for the Drude weight formula (|29j) with q ranging only over a = +,—■ 
As (|2*7jl has to be slightly modified by ± signs, because of e± = ±J^^-p' a , we find (|2*S)) 

consequently modified to p Q = ^27^ hi(l + r]' 1 ). The Drude weight is 

D= ii^Ly^ a r (ik^ni^) 2 fa (55) 

4tt/3 7 ^ /_ w (1 + ^(l + r/a 1 )^ In »7« 
Alternatively, in terms of a (= r/ + ), a (= r/Z 1 ) and due to the relation (JHU|) we get 

D =^Vf ( ^' n °° )2( *'7°' 2 dx. (50) 

with the notation a± := a, 02 := a. This expression with functions a and a calculated from 
(|39|) is our final analytic result for the finite temperature Drude weight within the spinon 
approach. In the next subsections we will numerically evaluate D for arbitrary temperatures 
and study analytically the zero temperature limit and the high temperature asymptotics. 



4.1 Numerical results 

Our numerical results are obtained for arbitrary temperatures T and various anisotropy pa- 
rameters A = cos 7 in the repulsive regime [0,1], cf. FigEJ We find qualitatively different 
behaviour depending on the value of A. 



For A close to the free fermion point A = the Drude weight D(T) is a monotonously 
decreasing function of temperature, see Fig|2 However, for larger values of A close to A = 1 
corresponding to the isotropic antiferromagnetic Heisenberg point the dependence of D(T) 
on temperature is non- monotonous! For sufficiently low T the function D(T) increases with 
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temperature. After taking a finite temperature maximum the function D(T) decreases. In 
particular, the values of D(T) at the isotropic point (A = 1) are non-zero for any value of 
temperature T, see Fig|3 This is in striking contrast to the results of jS] and those of the 
preceding section. A discussion of this will be given in the final section of this paper. 




4.2 Drude weight at T = 

In the limit of T — > and large argument x we have the scaling behaviour of the driving term 
in the NLIE 

^ ~n -||x| 



cosh | x 



(57) 



Prom this observation it was concluded jllj that the leading low T asymptotics of the physical 
properties is determined by values of the spectral parameter x ~ ±-ln/3. Indeed we find 
that limr^o a a (x ± — In 0) yields a well defined function of x still satisfying a non-linear 
integral equation. However the associated functions like ^a a and -§pa a satisfy linear integral 
equations. Actually, both sets of integral equations are identical up to different driving terms. 
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Figure 3: The Drude weight at the isotropic point A = 1 (7 = 0). A very steep slope close 
to T = is observed. This slope is infinite at precisely T = 0. 



Those, however, are strictly proportional to each other as they are given by 

d f f3 \ ^ T7r/3e -fH 
dx \cosh J 

8 



P 



dp \cosh 



(58) 



with proportionality factor ±§/3- Therefore, in the limit T — > and for the relevant range of 
spectral parameters x the derivatives of In a a with respect to x and P satisfy 



Hence in (|56JI we find the simplification 







J sin 7 



E 



In a a ) 2 (JL In a, 



(1 + a a )(l + a,, 1 ) 
By use of the "dressed function" formalism we obtain the identity 

/•oo q2 

12 ./o 4(1 + a a ) 2 



<ix = — /3 



dx. 



d 

[ln(l + a Q )]— e a dx, 



(59) 



(60) 



(61) 
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Hence the Drude weight is 

2 



( 7rsin7 
4 U~ 



J J E 7 Kl + «a)]eorfx, (62) 



where the term to the right of the partial derivative is identical to —2/3f. The second 
derivative with respect to j3h leads to the magnetic susceptibility 

1 vrsin7 
2u(7r — 7) 27 

Finally the zero temperature Drude weight is found to be 

Do = J. (64) 

871^ - 7) 

This result agrees exactly with the Drude weight directly obtained from the groundstate 
energy of a finite system with twist [U O2] • 

4.3 High temperature asymptotics 

As an illustration of the high temperature behaviour we show a plot of TD(T) for some value 
of A close to 1, see Fig0J The analytic evaluation using |)I8*|) in yields 

Da m, C (A)-^ (65) 

Again there is a problem with the analytically derived expression for C(A). Here it contradicts 
certain rigorously known properties of the high temperature limit that we are going to derive 
from the spectral representation (|2*|). First we symmetrize the summation 



-1 



I / „ \{n\J\m)\ 2 

> 



o.t f A - 



(66) 



^ 2^ : _^ K n l J l m )l 



2 £m 

> 

As D is the difference of two non-negative terms we see 

< ^<-r) (67) 

where equality holds if the second term in (|66[) disappears; this happens for A = where the 
matrix elements (n\J\m) for different states n and m yield zero. 

From direct calculations at high T we find (— T) ~ C /T with a constant C independent of 
A. From (|67|) we see C(A) < C/2L. On the other hand we know that equality holds in the 
case A = 0, therefore 

C(A) < C(0), (68) 

an inequality that has been obtained earlier on grounds of the optical sum rule, see e.g. p^. 
This is clearly violated by our analytic result l|65jl! 



The status of the high temperature results in the spinon approach and those in the TBA 
approach are equally problematic: in the spinon approach the inequality (|68j) is violated, in 
the TBA approach the symmetry with respect to A is violated. 
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Figure 4: The Drude weight at high temperature showing a 1/T-decrease (results shown for 
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5 Discussion 

In section 2 we reviewed the analytical method of computation of the finite temperature Drude 
weight as proposed in [7j and used in [S] for the study of the spin- 1/2 Heisenberg chain. We 
then employed this method as the starting point of our own analysis of the Heisenberg chain. 
This was done in two different ways. 

In section 3 we utilized the particle basis of magnons and their bound states. Within this 
approach we managed to reduce the resultant equations to only two non-linear integral equa- 
tions. These results are totally equivalent to those of [Hj, however, they allowed for an analytic 
continuation to the treatment of all anisotropies — 1 < A < 1. Also, we derived an analytic 
formula for the high temperature asymptotics. 



In section 4 we employed the spinon and anti-spinon particle basis. The results obtained in 
this approach strongly deviate from those of the preceding section and hence from [S] . We ob- 
served qualitatively different behaviour of D(T) for A close to and 1 showing monotonous 
and non-monotonous temperature dependence, respectively. Instead of a drop of D(T) at 
T = our results show an increase with apparently infinite slope at T = 0, see FigJHl This 
is reminiscent of the behaviour of the magnetic susceptibility with infinite slope due to loga- 
rithmic corrections at T = 0. 
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Though disagreeing with [Sj we find strong similarities of our results in the spinon approach 
with the numerical work ^S] which is based on complete diagonalization of quantum chains 
up to length L = 14. Unfortunately, the non-monotonous dependence of the Drude weight 
data obtained in |16j may still be considered as plagued by finite size corrections. In a later 
numerical analysis in JT£| the size dependence was carefully studied for system sizes up to 
N = 18. It was argued that a non-monotonous behaviour found at A = 1 was not an 
artefact of finite size effects and qualitatively agreed with our (then unpublished) data that 
we presented here in section 4. Also, by use of quantum Monte Carlo calculations at low 
temperatures, the authors of ^H] found quantitative agreement of their results for anisotropy 
A = cos 7r/6 with our (then unpublished) findings in the spinon approach, which are in strong 
disagreement with jS]. 

Still, the numerical treatments have to be considered with care as the size dependence is 
indeed very strong. However, also the analytical treatments have to be considered with care. 
For the case of the TBA approach on the basis of bound states (strings) we have indicated 
certain problems and mechanisms of possible failure. In our understanding, signatures of this 
are even visible at high temperatures, cf. subsections 3.2 and 4.3. This will be studied in 
more detail in |13j . 

In the case of the spinon approach we are confident that the results are reliable at low temper- 
ature, whereas they are not trustworthy at high temperatures. The reason for this is simply 
that the concept of the spinon and anti-spinon particles is field theoretic and restricted to 
low energies and low temperatures. (Strangely, for the static properties encoded in the free 
energy, the spinon-concept gives correct results for all temperatures and fields! Some further 
aspects along this line are discussed in the appendix.) 

Finally, we like to mention that our results in the spinon approach have been confirmed 
by conformal field theoretical arguments developed in ^5] for certain anisotropy values A, 
however for some other values (notably those close to 1) there is still strong disagreement. A 
comparison of all data in collaboration with the authors of ^5] is on the way. 
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Appendix 

Here we address the question whether the relations (|50[) and (|53jl are more than accidental. 
We want to argue that the functions p a (a = ±) obtained from equations ()28|) in the spinon 
approach (and satisfying (j9)lljl ) are true density functions. A (necessary) criterion to be 
satisfied is the requirement that p a determine all thermal expectation values of the higher 
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conserved quantities F n where 



F n = |^logT(x) , (A.l) 

OX n x=0 



and T(x) is the six-vertex model row-to-row transfer matrix in dependence on the spectral 
parameter x with decoupling point x = 0. For n = 0, 1 we have Fq ~ P (momentum operator) 
and F\ ~ H (Hamiltonian). We are interested in the cases n > 1 and define the differential 
operator D n = d n /dx n . 

The contribution of a magnon with spectral parameter x to the eigenvalue of any F n is 
(D n p)(x). If p a determine all thermal measurements we must find the relation 

/oo 
(D nPa )(x) Pa (x) dx. (A.2) 

We want to verify this expectation and calculate 

(F n ) = (D n logT(x)) = ^ logTr [exp(l + zD n logTje"^] , (A.3) 

where the derivative is to be taken at z = 0. The reason for introducing the rather involved 
expression on the r.h.s. is that exactly this quantity can be calculated in the QTM approach 
resulting into almost literally (|54[) except the replacement 

0e a (x) -» 09 - zD n „ x )e a (x), (A.4) 

which can be derived by simple but lengthy calculations. Note the subscript n — 1! The free 
energy (or more precisely the quantity logTr[...] in (|A.3|0 is still given by the r.h.s. of (|53|) . 
Hence the thermal average of F n is 

The expressions l)A.2|) and (|A.5|) look quite different. Still they are equivalent as can be seen 
in the "dressed functions" formalism. To this end we note (|54jl with the replacement (|A.4jl . 
and take the derivative with respect to z at z = 

d Id 
(1 + r) a )— In (1 + t?^ 1 ) = -D n ^e a + — ^ x<*p * tt In (l + J?^ -1 ) • (A.6) 

Z IT ^ Z 

This equation can be regarded as an integral equation for ^^(l + r/^" 1 ) and is very similar 
to that one for p a Q, just the inhomogeneity (driving term) is different 

(1 + Va)pa{x) = ^-Pa( x ) + E Xa P * Pi 3 ' 

P 

It is now a standard exercise in mathematical analysis to show the identity of (|A.2|) and (|A.5|) . 



Unfortunately, the criterion that the functions p a yield the correct expectation values is 
necessary, but not sufficient. True density functions are usually real and non-negative. In our 
case, however, complex valued functions are in principle allowed as we may have distributions 
along curved lines in the complex plane. Along these lines just the ratio of particle and hole 
densities (i.e. the function rj a ) must take positive values. In the spinon approach this is 
satisfied for low temperatures, but not for high temperatures. 
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